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In this paper a new representation formula for the spherical components of displacement 
in the theory of elasticity is obtained. This formula may be used to reduce mixed boundary 
value problems for the elastic sphere to standard problems in potential theory. 



1. Introduction 

In treating boundary value problems in the theory of elasticity it is frequently convenient 
to represent the components of the displacement vector in terms of harmonic functions, and 
thus to reduce the elasticity problem to one in potential theory. Several such decompositions 
in terms of harmonic functions are known, [7, p. 342-343], 3 [1, p. 231], [4], the most widely used, 
perhaps, being the Papkovitch-Neuber decomposition (see for instance [7, p. 331]). 

Half-space problems have been quite thoroughly treated from this point of view, as in 
fact have the first and second boundary value problems for the sphere. However, in treating 
mixed boundary value problems for the sphere, the standard decompositions are not suitable. 
Such mixed problems could be simplified considerably if it were possible to obtain convenient 
representations for the spherical components of displacement in terms of harmonic functions. 
We obtain such a representation in this paper and demonstrate how it yields in a simple man- 
ner solutions to standard mixed boundary value problems for the sphere. 

In this paper whenever an index is repeated in a single term, summation from 1 to 3 will 
be implied. Use of the symbol /,< to denote c)//dx z -, where x t is one of the rectangular coordinates 
(#i, x 2 , x 3 ), will be made throughout. 

Let D be a three dimensional region which contains an open subset Q of a spherical surface 
of radius a with center at the origin of a spherical coordinate system (r, d, cp). The region D is 
assumed to be a portion of an isotropic homogeneous elastic medium. 

We suppose that D does not contain the origin. Furthermore, we assume that every ray 
from the origin which intersects D, intersects it in a single line segment. 

In D the component of the displacement vector satisfies the equilibrium equations 

Ui tjJ + aUj tji =0 (1.1) 

where a involves the elastic constants (a=(l — 2a)~ 1 where a is Poisson's ratio). As far as our 
representation formulas are concerned we need put no restrictions on the constant a other than 
that it be different from —1. However, all physically interesting values of a lie in the interval 
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In this paper we obtain the following representation for the spherical components of dis- 
placement in D: 

x iUj =r> ^2G+(2-a)^G/dr-ar^( K r^y+r^ (1.2) 

tijkXjW k =r 2 -I —(2+3a)e i j k XjG, k —ar^-e ijk x j G, k —(H+r-^ 7 J V 

where x//, G, and Hare harmonic functions in D. Consequently all terms in braces in (1.2) and 
(1.3) are also harmonic. In (1.3) e ijk is the permutation symbol 



*ijk- 



if any two indices are equal 

1 if i jk is an even permutation of 1 2 3 (1-4) 
— 1 if ijjc is an odd permutation of 1 2 3. 



2. Representation Formulas 
We first rewrite (1.1) as 

€lJk dx 6 ^m^l,m+(a+l) U jtJl = 0. (2.1) 

For a 9^ — 1 we'see from (1.1) that Uj,j is harmonic. It has been shown that the representation 

%,,=(l+^)(6ff+4rg), (2.2) 

where G is harmonic, is valid in D (see for instance Bramble [2]). Equation (2.1) then may be 
integrated to give 

e<i*tt/.*=(a+l) *ijkXj(6G+4r^y+x,i (2.3) 

where x is a harmonic function. We express x hi D as 

X=$+r^- r - (2.4) 

Thus Uj is the component of a vector whose divergence is represented by (2.2) and whose curl 
is expressible by (2.3). The component Uj may clearly be expressed as the sum of a particular 
solution of (2.2) and (2.3) and the gradient of a harmonic function, i.e., 

u i =2x i ^G+(a+2)r^'j-r^2(l+2a)G, j +ar^ r G.^+e itk X^, k +t i (2-5) 

where \p is harmonic in D. 

Using (2.5) we now form the expressions for XjU 3 and € ijk XjU k which are respectively the 
scalar product and the vector product of the radius vector and the displacement vector. 

^"•O+e-Org-r |(rg)}+r£ (2.6) 

{d 1 d$ 

— (2+3:*) e Uk 3: j G tk —ar^ t € ijk x i G 9 j c —<l> i f + W-fo +«<j»*/A*- ( 2 - 7 ) 
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and 



Equation (2.6) is precisely the desired expression (1.2). In order to obtain (1.3) we substitute 
for<f> in D (see Bramble [2]), 

$=tf+2r^ (2.8) 

where H is harmonic. We then write 



dr 



+{-£+*!('£)}-•■('£)>('£)... 



As noted previously, the term in square brackets is harmonic. Equations (2.7), (2.8) and (2.9) 
yield immediately the desired representation (1.3). 

In obtaining our decompositions we have frequently made use of the fact that in the 

region in question any harmonic function could be represented as ( 1 -\-d\r ^— ) ( l-\-a 2 r >r ) 9 

where g is harmonic. We have used only representations in which a\ and a 2 were non-negative. 
Bramble [2] has shown that such a representation is valid in D regardless of the sign of a x and 
a 2 . However, for a x and a 2 non-negative it is easily verified that such a representation is also 
valid in any region whose boundery is star-shaped with respect to the origin. Consequently 
(1.2) and (1.3) are valid in any such star shaped region. In the next section we use the repre- 
sentations (1.2) and (1.3) to solve certain boundary value problems for the elastic sphere. 

If the elastic region is exterior to a surface C which is star-shaped with respect to the 
origin and if the components of the displacement vector are 0(r~ 2 ) as r-^co } it can again be 
shown that the decompositions (1.2) and (1.3) are valid therein. 

3. Mixed Boundary Value Problems for the Sphere 

Let u t satisfy (1.1) on the interior of an elastic sphere. On the surface r=a let 

y.- l a x T i jX i x j +a 2 x j u j =f(e } (p) 

(3.1) 
m -1 &i UjkTkiXjXi+b 2 € ijk x j u k =g i (e ) <p) 

where T tj is the stress tensor, /x is the shear modulus, and where/ and g t are prescribed boundary 
data. Equations (3.1) may be written equivalently as 

ai(a+l)/-^Ofy^) + [ai(a— ^)+a 2 ]x J Uj+a 1 (a— l)e iJfc x i ^-€ klm XiU m =f(e j <p) 

bir^e ijk x j u k +(b 2 —2b 1 )e ijk x j u k +bie uk x J (x l u l ) tk =g i (d,<p) (3.2) 

on r=a. Using (1.2) and (1.3) we can actually obtain the following representation for the 
left hand side of (3.1) in terms of harmonic functions 



jtx Wijtt 



dr 



i x j +a 2 XjU j =r 2 < — 2aia(rr— j G—a 2 a(r^-) 6r+[2ai(4a:— l)+a 2 (2— a)] r ~^— 

+2[(3a-l)ai+a 8 ]fl f \+2flir ^ (r ^j-{2a l -a 2 )r 
^ l bie i j k r kl x j x l + b 2 e ijk x J u k =rh ijk Xj^^—2ab 1 (r G— (6 2 + 4&i)ar^ + [26i--(2 + 3a)6 2 ]G f i 

_ 2r ,( 6>r |. + i,)(H +f ^) ( +,„ A A{2J I ^ + (J,-26,)*.} 

(vA +i ,- 2S ,){4,.^ + ^(,^)]_,(»)J. ( , 3) 
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From (3.3) we obtain expressions of the form 

-<•' [(' +' I) (*" ^+«0].,= P W 

valid throughout ZX The constants A t and B t are determined immediately from (3.3) and 
P(j) and P(gi) denote Poisson integrals. 

Forming XiP{g^) we obtain the expression 

Under suitable conditions on B & and B 7 this expression may be integrated to give an expression 
for H. With H known throughout D, equations (3.4) give two identities in the two unknowns 
G and \p. The function \p can be eliminated by taking the proper combination of r derivatives 
of the first of equations (3.4) and r derivatives of the expression formed by operating on the 

second of equations (3.4) with € Zmz ^^ — This results in an equation involving a linear 

ox m 

combination of r derivatives of G. If this equation can be solved for G then the insertion of 

G into the first of equations (3.4) yields the expression for \p which completely determines the 

displacement vector in D. It is a simple matter to obtain a formal solution for G, for we have 

an equation of the form 

§o,(r|)W. (3.6) 

We expand the known harmonic function F in terms of spherical harmonics, i.e., 

2 ^ (r j£) G=ib S r n P% (cos 6)[A nm cos m<t>+B nm sin m# (3.7) 

i=o \ or/ w =o m=o 

The solution to this equation is 

oo n [~ N "1-1 

G=J2 S S ^ r w P™(cos 0) [Am cos m4>+B nm sin ra</>] (3.8) 

71 = 771 = L*' = o J 

provided 

S^VO (3.9) 



i=0 



for w = 0, 1,2, . . ., or any positive integer. For a given a, (3.9) imposes certain restrictions 
on the constants a h a 2 , &i, and b 2 . 

As a particular example consider the case in which a 2 and b } in (3.1) are zero. The constants 
#i and b 2 can without loss be taken as unity. Thus equation (3.4) becomes 

^[-2a(r^)V+2(4a+l)r^+2(3a-l)(?]+2r| ; (r|;-l)^=P(/) (3.10) 
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and 

aUi *&* dx k l ar £^"~ (2+3«)6 ? J- e iJk Xj^ k -a 2 [H+r^rJ % 

where P(g-D and P(f) denote Poisson's integrals. If we take the scalar product of the radius 
vector and (3.11) we obtain the differential equation for H 

which yields the result 

H ~J P {p 2_ a 2 ) d P . (3.13) 

It is possible to add an arbitrary constant to the right hand side of (3.13) but since only deriva- 
tives of // occur in (1.3) we may take the constant equal to zero. We now operate on (3.11) 

with the operator e lmi Xi r- — to obtain the expression 

r 5P {r[*-(2+3«)«*7-awg] } = a,,A^-P(y t ) (3.14) 

which yields for \p 

^a^+^G+ar^-^^p-^r-tiu^—PigMp- (3.15) 

Again we have set the arbitrary additive constant equal to zero in (3.15) since only derivatives 
of \p occur in (1.2) and (1.3). Substitution of (3.15) into (3.10) yields the result 

wlicre 

Cl ={(3-a) + [17-22a-23a 2 ] 1/2 }[4(a+l)]-' 

(3.17) 
C2 ={(3-a)-[17-22a-23a 2 ] 1/2 }[4( a +l)]- 1 

and 

R=P(f)+2 (r |.-l) [r-'J^mXi ^ P(g*)d P ] (3.18) 

For a>l/3 (which corresponds to the range of a which is of physical interest, i.e., — l<V<0/2) 
C\ and c 2 have real part less than unity. Tims if we operate on (3.16) by rd/dr we obtain 

('i-*)('i-«>f-t««+ 1 >i"--f <"» 

an expression which can easily be integrated for G, giving 

G=[4a>(«+\)}- 1 { g( y } + J V'- 1 J] V"- 1 £ Pi"' ^ *wMp}" (3.20) 
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If the 


indicated integration h 


1 performed we 


obtain for 


Cit^c 2 the 


expression 








G=l4a*(a+1)}- 


(^ C1C2 


-|H> 




-sr>- 


,55 
dp 


*]}■ 


(3.21) 


If c,= 


-c 2 (which means 


that 


a satisfies the 


equation 


17— 22a- 


-23a 2 =0) 


(3.20) yields 


upon 


integration 



















GH^cKa+l)]- 1 |<vr c 'JY Cl log {rip) || ^P-^ P~ Cl |f *>+^}' ( 3 - 22 ) 

The function ^ is now determined by insertion of (3.21) or (3.22) in (3.15), and the problem is 
solved. 

4. Concluding Remarks 

In a recent paper [3], the authors indicated how to obtain an explicit continuation across 
a spherical surface Q on which/ and g t of (3.1) were zero. The decomposition (1.2), (1.3) could 
also be used to obtain this continuation formula. It is also possible to obtain in a straight- 
forward way from (1.2) and (1.3) sphere theorems with various combinations of surface trac- 
tions and displacements vanishing on the surface of the sphere. 

If in (1.2) and (1.3) a is allowed to tend to infinity while at the same time aG tends to a 
harmonic function F, then these decompositions may be interpreted as representation formulas 
for the spherical components of the velocity vector in the slow motion of viscous fluids. (See 
for instance [6, p. 522].) 
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